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^ ■ Abstract 

A recently constructed strangeness-including Statistical Bootstrap Model (S-SBM), 

D , which defines the limits of the hadronic phase and provides for a phase beyond, is fur- 



ther extended so as to include a factor 7^ that describes strangeness suppression. The 



^ . model is then used to analyse the multiplicity data from collision experiments in which 

d ' the colliding entities form isospin symmetric systems, the primary focus being on 5-|- 5 

interactions (NA35 collaboration). An optimal set of thermodynamical variables is ex- 
tracted through a fit to both the inclusive 47r and midrapidity data. The assumption 
that the measured particles originate from a thermally and partial-chemically equi- 
librated source described by the S-SBM is satisfactorily established. The proximity 
of the thermodynamical variables extracted from the S + S data to the limits of the 
hadronic phase is systematically investigated. Finally, experimental data from p + p 
collisions (UA5 collaboration) are similarly analysed. 



PACS numbers: 12.40.Ee, 05.70.Ce, 05.70.Fh, 25.75.Dw 



1. Introduction 

The description of multiparticle hadronic states which emerge from a coUision process 
at very high energies can, most economically, be conducted by adopting thermalization 
hypotheses pervading their evolution course. It is, in fact, the quintessence of the statistical 
mechanical approach to a given system with a large number of degrees of freedom that it 
transcribes a hopelessly complicated microscopic description into a corresponding concise one 
which employs a handful of thermodynamical variables (e.g. temperature, volume, chemical 
potentials and the like). 

Naturally, the ultimate goal of any attempt to address issues involving multiparticle 
production in high energy collisions is to account for the dynamics which give rise to the 
experimentally recorded profile of the final state of the system. For the particular case of 
heavy ion collisions, on which the bulk of our attention will be focused in this work, the 
question of fundamental interest is whether the overall collision process has gone through a 
deconfined phase and, if yes, in what way can one pick signatures of this occurrence in the 
composition of the final state. 

Given the above general remarks, let us now become more concrete and consider the basic 
thermodynamical scenario which has been adopted by the pioneers of this approach [1-3]. To 
describe a system that has been produced by one (or more) high energy density source(s), one 
assumes a dynamical evolution that reaches freeze-out points (thermal, chemical) at which 
corresponding thermodynamical quantities acquire their final equilibrium values, reflected in 
the observed particle species and multiplicities. 

The statistical mechanical analysis evaluating the experimental data can be conducted 
in the spirit: (a) Taking for granted the validity of the thermal picture and looking for 
possible discrepancies which provide signals of something interesting having taken place, (b) 
testing the reliability of a given thermal model, (c) a combination thereof. In any case, it is 
particle multiplicities (or ratios) which provide concrete numerical input for the theoretical 
analysis. Such information addresses itself to the chemical composition of the system, which 
means that chemical potentials become a key ingredient in this type of description. In 
this context, one must also consider the prospect of amending the equilibrium scheme by 
inserting parameters which account for partial equilibrium conditions with respect to a given 
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quantity. This possibility will, in fact, prominently enter the present analysis in connection 
with saturation of strangeness phase space. 

Theoretical confrontations of thermal hadron production in various types of experiments, 
ranging from e^+e~ to nucleus-nucleus collisions, have been conducted within the framework 
of a statistical model corresponding to an ideal hadron gas with results that, a posteriori 
at least, justify the premises of the whole approach [1-11]. However, such a model has 
a built-in assumption that the hadronic world has no bounds and it does not anticipate 
the existence of a new phase, beyond the hadronic one, when a critical point in the space 
of thermodynamical parameters is reached. By contrast, the statistical bootstrap model 
(SBM), originally introduced by Hagedorn [12], see also [13-16], does anticipate the end of 
the hadronic domain and, under certain conditions on model parameters [17], is consistent 
with the existence of a new phase beyond the hadronic one. Given this state of affairs, as 
well as the fact that the quantum number of strangeness plays a central role in revealing 
the quark-gluon plasma (QGP) phase, we have recently constructed an extension of the 
SBM so as to incorporate Strangeness [18,19], henceforth to be referred to as S-SBM. The 
model exhibits attractive properties consistent with what one anticipates from the transition 
between the hadronic and the new, presumably QGP, phase of matter. 

At present, the S-SBM treats u and d quarks on the same footing so that a single fugacity 
corresponds to both of them. This amounts to considering only isospin symmetric systems, 
i.e. systems for which the number of u-quarks minus u-quarks is held equal to the number 
of d-quarks minus d-quarks. Then it is easy to prove that 

riu - riu = rid - nj ^ = = Xq , (1) 

where the n's stand for quark numbers and the A's for corresponding fugacitie^. The above 
relation, in the presence of the condition that the total strangeness vanishes, is equivalent to 

B = 2Q^Xq = 1 , (2) 

where Q is the total electric charge and Xq is the corresponding fugacity. So in the afore- 
mentioned situation we need not consider the existence of the additional fugacity Xq and 
the present form of S-SBM ascribes precisely to these systems. 

^Our notational conventions regarding fugacity labels will be specified in the next section. 
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In this paper we intend to confront, primarily, experimental data coming from isospin 
symmetric S + S collisions (NA35 experiment at CERN) within the framework of the S- 
SBM. For purposes of completeness, we shall also use our methodology for the evaluation 
of p + p data (UA5 experiment at CERN). We shall be addressing ourselves to both hadron 
multiplicities and corresponding particle ratios which involve strange particles, assuming 
partial strangeness saturation. To this end, we shall introduce into our scheme the additional 
parameter 7^, which will be disscused in the following section and can be identified as the 
"fugacity" pertaining to the number of strange plus antistrange quarks. 

Our analysis will cover both the inclusive Att set of data and the, more restricted, midra- 
pidity region. A fine tuning due to corrections from Bose/Fermi statistics will also be taken 
into account, given that the SBM scheme adopts the Boltzmann distribution (maximum en- 
tropy content). Our purpose is twofold: first to find out whether a thermal description based 
on the S-SBM accounts for the experimental data in a satisfactory manner, and second to 
find the proximity of the data to the critical surface which sets the limits of the hadronic 
phase. 

The organisation of the paper is the following. In the next section we shall consider 
the ramifications brought by the inclusion of the 7^ parameter into the S-SBM scheme. In 
Section 3 we shall proceed to present the methodology by which we propose to conduct our 
analysis of observed particle multiplicities. The actual examination of the experimental data, 
both for the S + S and p+p cases, will be presented in Section 4. A fit will be employed to 
asses our theoretical predictions. We shall devote the last section to a discussion-evaluation 
of our results, paying special attention to the issue of whether and/or how close the S + S 
data come to revealing the attainment of the QGP phase, always according to our model. 
Two appendices (A and B) are devoted to technical matters concerning routines used in 
our numerical calculations, while in appendix C we discuss the particle ratios vs particle 
multiplicities issue. 
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2. Partial strangeness equilibrium in the S-SBM 

We shall begin this section by presenting, in an outline form, the S-SBM construction of 
Ref. [18,19]. The underlying assumption in that work was that the dynamical development 
of the hadronic system formed after a high energy collision is characterised by chemical 
equilibrium throughout. Subsequently, we shall consider the generalisation of this scheme 
brought about by allowing for only partial strangeness equilibrium. 

The two basic ingredients going into any SBM scheme are the following: 

(a) The so-called bootstrap equation, which incorporates in a self-consistent manner the 
internal dynamics operating in a relativistic multiparticle system. 

(b) A statistical mechanical account of the said system in terms of a suitably defined 
partition function (equilibrium mode of description) 

The bootstrap equation refers to a spectrum function and receives input from the set of 
all known hadrons through the so-called input function. Relativistic invariancc in the casting 
of the system leads to a term of the form B{p'^)t{p'^, b, s), where B{p'^) is a kinetic factor and 
t(p^, 6, s) the naturally defined spectrum function which exhibits dependence on baryon (6) 
and strangeness (s) numbers. Rearrangements of the kind 

B{p')T{p\b,s)^B{p')f{p\b,s) (3) 

define different versions of the bootstrap scheme depending on whether the f-function carries 
all or part of the dynamics acting internally. In the latter case, the ^-factor can also assume 
a dynamical role. 

The particular choice specifying the S-SBM has been extensively discussed in [18,19], on 
the basis of the decisive physical advantages it exhibits [20] , is 

(27r) (27r) 

where the last expression refers to the rest frame of the particle/fireball. The above selection 
simply implies that B(jp), which takes the form B{m?) in the rest frame, remains a purely 
kinematical quantity. 

Setting B{m^) = Hqui^ one arrives at the notable relation 



which hnks the bootstrap scheme parameter Hq with the MIT bag constant B, the latter 
entering through the relation = It turns out that Hq is directly proportional to the 
maximum value of the critical temperature Tq. 

A Laplace transformation which takes the set of variables (p^, b, s) to the set Xb, Xs) 
introduces the inverse four-temperature, along with baryon and strangeness fugacities, into 
the bootstrap scheme. One is thereby led to the construction of a partition function 
Z{P,V, Xb, \s) in the fireball rest frame. In this way a statistical mechanical description 
of the system is introduced, giving the final form to the S-SBM. 

The analysis conducted in Ref. [18,19], with respect to the S-SBM, took place in the 
3-dimensional space of the thermodynamical variables {T, /ig, /ig), where /ig and /ig stand for 
up-down and strange quark chemical potentials. In this space, the following two surfaces 
were considered: 

(a) The critical surface specified by one of the conditions 

^{Tcr, fJ'qcr, fJ'scr't Hq) — la4: — 1, G{Tcr, fJ>q cr, fJ's cr'i Hq) — la2 , (6) 

where (p corresponds to the input and G to the mass spectrum-containing function entering 
the bootstrap equation. This surface constitutes the earmark of the bootstrap scheme, which 
separates it from the ideal hadron gas model. It signifies the termination of the hadron phase, 
since the bootstrap equation does not posses a physically meaningful solution beyond the 
critical surface. In other words, the critical surface sets the boundaries of the hadronic world 
in the space {T, fig, fig). 

(b) The < S >— surface which imposes a strangeness neutrality assumption on the S-SBM 
construction and is specified by 

p^aG(...KX,H„)^^ = . (7) 

J/3 oXs 

As already mentioned in the introduction, our aim in this paper is to extend the frame- 
work of the S-SBM so as to allow for partial strangeness equihbrium before confronting 
experimental data (particle multiplicities or ratios). This is most conveniently done by in- 
troducing the variable 7^ which is, in fact, a fugacity related to the number of s-quarks plus 
s-quarks [21] (we shall henceforth call this number |s| = rig + ris): 

7s = X\s\ = exp(/i|,|/T) , (8) 
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The bootstrap equation in this generahsed model reads {gbs\s\ denotes degeneracy factors 
apphcable to the given set of labels) 

oo -1 „ n 

BipW, b, s, \s\) = 9,s\s\BipW - ml\s\) + - EpO- 

V n=l 1=1 

input term 

TL fL TL Tl 

■E^K{b-^k)Y: Sk{s - E «0 E ^k{\s\ - E klO n B{pMpI K Si, \s\i)d% . (9) 

{bi} i=l {s,} i=l {\s\i} i=l i=l 

In the above equation, baryon number is denoted by "6" and strangeness by "s" . Then we 
can perform in (9) four Laplace transformations which lead to the replacement (after going 
to the rest frame of the system) 

{jf,b,s,\s\)^{P,XB,Xs,%) (10) 

where Xb and Xs are the fugacities corresponding to baryon number and strangeness, re- 
spectively. Since we are accustomed to working with fugacities corresponding to quarks, we 
can equivalently use the set defined by 

x,^xT, xr^isxTxs\ Ar=w/%. (11) 

The up and down quark fugacity is denoted by A^, whereas A^^* and Aj*^^ denote the s- 
quark and s-quark fugacites respectively. The index neq means that the chemical equilibrium 
connected with strangeness is, in principle, not acheived. The factor 7^ entering (11) can, 
then, be given by 

In order to stay in agrement with the conventions of the quark fugacities used in [18,19], 
as well with the ones used by other authors, we shall employ the set (Ag,As,7s). Now 
Xs = A^^^A^^ is the fugacity of the s-quark which corresponds to the chemical equilibrium 
of strangeness. 

The bootstrap equation acquires the form 

(^(T, Ag, Xs, 7j = 2G(T, Ag, A^, 7J - exp[G(T, A^, A^, 7^)] + 1 , (13) 

with the functions (p and G given by 

<f{T, Xg, Xs, -fs; Ho) = 27rHoTY,KT.9siml,Ki{m,i/T) , (14) 



G{T, Xg, A„ 7,; Hq) = 2ttHoT / mVoK, A^, A„ -fs)K,{m/T)drrv' , (15) 

Jo 

where the subscript "o" on r signifies the S-SBM choice given by eq. (4). The fugacity Aa 
corresponds to the existing hadronic famihes that are used as input in (14) (fight unflavoured 
mesons, kaons, A?" & A, A & E, S and Q Baryons). It runs over all particles and antiparticles 
and obeys the equality 

where rii corresponds to the number of the ith flavour quark which is included in the "a" 
hadron. For the case of the light unflavoured mesons which have the quark content of the 
form ciqq + C2SS (ci + C2 = 1) the corresponding fugacity equals 

Aa = Ci + C275^ . (17) 

For the evaluation of the coefficients ci and C2 we have used [22]. At the same time the 
partition function is amended by the addition of the fugacity variable A|s| or, equivalently, 
the chemical potential ii\s\. 

The critical and < 5" >= surfaces are now determined by equations similar to those 
given by (6). The only difference is that 7^ is also included as a variable: 

'f{Tcr, fJ'q cr, fJ's cr. Is cr'i Hq) = In 4 — 1, G{Tcr, Hq cr, fJ's cr. Is cr] Hq) = In 2 , (18) 

and 

r,.8G(.,X X.,j,,H„)^^ = . (19) 

J(3 OAs 

For the S-SBM, as constructed in Refs [18,19], Tq, the critical temperature for vanishing 
chemical potentials, can be taken as the only free thermodynamical parameter of the model, 
in the place of Hq. Tq is determined by the relation 

(/?(To,A, = l,A, = l;//o)=ln4-l, (20) 

and, as can be seen, is function of Hq only. 

In the present case the last equation becomes 

(^(To, A, = 1, A, = 1, 7,; //q) = ln4 - 1 , (21) 
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and To is a function not only of i^o, but of 7s as well. It should be noted, on the other hand, 
that the bag constant B continues to be in one to one correspondence with Hq, as can be 
seen from (5). 

The dependence of Tq on 7^ complicates the numerical analyses needed to determine the 
profile of the critical surface in the considered extension of the S-SBM. Nevertheless, the 
relevant study has been carried out and the results are presented in a series of figures. 

Figure la depicts two characteristic intersections of the critical surface. One is with the 
Us = MeV and the other is with the fis = 500 MeV plane, for different values of 7^. The 
unphysical value 7^ = 1.50 has been chosen to show the effect of the increase of 7^ beyond 
unity. Figure lb depicts the same situation as the previous one, but the intersections of 
the critical surface are with the planes /ig = MeV and fig = 300 MeV. As it can be 
seen, a decrease in the value of 7^ below 1 leads to an expansion of the critical surface. 
The opposite happens for the increase of 7^ above 1. Figure 2 shows the variation of the 
critical temperature Tq with 7^. The decrease of 7^ leads to the increase of Tq and vice versa. 
Figure 3 displays the < S >= surface for the S-SBM and for the "ideal hadron gas" . The 
corresponding merging with the critical surface, in the S-SBM case, for two characteristic 
values of 7^ (1 and 0.5) is also depicted. The above curves are plotted for constant fugacity 

(/i, = 0.4). 

In figures 1-3 a specific value of Ho has been used which leads to Tq = 183 MeV for 7^ = 1, 
corresponding to the maximum value of the bag constant, B^^^ = 235 MeV, according to 
[19]. By comparison, figures 4a and 4b present the intersections of the < S >= with the 
critical surface for two values of Hq (one leads to Tq = 150 MeV and the other to Tq = 183 
MeV, for 7^ = 1). Fig. 4a shows the projections of the aforementioned intersections on the 
A^s ~ f^q plane, while Fig. 4b the projections on the fig —T plane. 

In conclusion, the extension of the S-SBM to allow for partial strangeness equilibration is 
feasible and its impact on the model is simply to shift the critical and < S >= surfaces as 
7s changes values. Its inclusion in the S-SBM scheme provides an extra (thermodynamical) 
degree of freedom, widening its scope of applicability. 

^The strangeness production is usually found to be suppressed, so physical values for 7^ should be less 
than 1. 
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3. Methodology for Examining Experimental Data 

Our confrontation of experimental data has as its focal point particle multiplicities. Par- 
ticle ratios will also be studied at a subsequent stage. Our task will be to determine the 
thermodynamic parameters (V^, T, A^, A^, 7^) which best fit the {N) particle multiplicities 
measured in a given experiment. In principle, this is equivalent to determining the parame- 
ters (T, Aq, As, 7s) which best fit a group of (A^ — 1) indepedent particle ratios we can form 
from the N measured multiplicities. Although more complicated, we shall prefer the first 
case, where it is possible, for the reasons discussed in Appendix C. 

Introducing hadron fugacities Aj into our partition function via an extension of the form 

Z{V,(3,Xg,K,'ys) ^ Z(V-,/5,A„As,7a,...,A„...) , (22) 

we have 



j^thermal ^ ^ d\nZ{V,P , \, As, 7s, . . . , A^, . . 



(23) 



or equivalently 



A^f~'(^,T,A,,As,7s) 

VT'^ i-T 1 1 (993(y,Ag,As,7s,...,Ai, ...) 



dy(24) 



=1 



A-k^HqJq |/5 2-exp[G'(7/,Ag,As,7s)] dXi 

where TV/'^'^''"*"' stands for the number of particles of type i coming directly from the collision. 
To this we must add the number of particles of the same species resulting through secondary 
production processes. Accordingly, the number tv*'^^'"^^, representing the theoretical predic- 
tion, is given by 

■^theory jythermal _j_ thermal (25) 



where the bij are are branching ratios corresponding to the decay of resonance j into a 
particle of type i. 

As can be seen from (24), the system's volume enters in each particle multiplicity as a 
common multiplicative factor. This does not change if finite volume corrections are intro- 
duced to account for particle size, as in [14] . In such a case, the volume V entering in (24) 
should be replaced by the free volume A, given by 

V 

^ ^ l + e,,iT,{X})/4B ■ ^^^^ 
9 



In the above relation, Spt is the energy density of a system of point particles which is a 
function of the same variables (T, {A}) as those characterising the system with the extended 
particles ({A} denotes all fugacities collectively). 

The actual determination of V itself is of no importance to our analysis. It will only 
be viewed as a parameter which serves to fit the data. So we need not consider volume 
corrections. Moreover, the analysis can be simplified (e.g. evaluation of the first and second 
derivatives with respect to the temperature) if we choose VT^ /Att^ as our free parameter in 
the place of V . We shall, therefore, adopt this choice in the following. 

Our first task is to determine an optimal set of values for the thermal parameters, i.e. a 
collection {VT^ /An'^, T, \g, Xg, 7^) which best fits the experimentally measured multiplicities. 
A primary concern is the imposition of the constraint < S >= 0, which, according to the 
S-SBM scheme, amounts to enforcing eq. (7). This constraint does not allow, of course, the 
consideration of all the above parameters as being free (one of them has to be determined by 
the rest). Our way of ensuring the constraint is by introducing an extra parameter / which 
plays the role of a Lagrange multiplier. 

Putting the above considerations together we form the function^ 



^iVTV.A T, A,. A.. ... ^ E rivr - iV--(^T3/4.3, r, A„ A„ 



1 2 



i=l 

^ir^s dG{x,X,,Xs,ls) ^^ (27) 

J/3 OAs 

where i runs over all hadrons measured in the experiment and cij is the corresponding experi- 
mental error. Our aim will be to determine the values of the parameters 
{VT^/Atc^, T, Xq, Xs, 7s, /) which minimise the above function. The problem amounts to solv- 
ing a system of six equations of the form 

^^%-^ = (, = 1.....6). (28) 

where the Xi run over the set of parameters (VT'^/iTi'^, T, Xq, Xs, 7s, /) 

The above system can be solved with a generalisation of the Newton-Raphson method 
to a multidimensional space. By this procedure we have to evaluate the second partial 



^We have denoted the function by for obvious reasons: It provides the measure of a chi-square fit. 
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derivatives of with respect to its parameters. The method is quite sufficient when the 
point, which represents the optimised parameter values, is well inside the domain of the 
hadronic phase. On the other hand, it becomes very ineffective when this point is near 
the critical surface (or, worse, outside)^ This occurs because, when we begin to search 
for the solution with the Newton- Raphson method by giving an initial starting point, the 
subsequent points, through which the function passes during the evaluation procedure, in 
general oscillate. In the multidimensional space this oscillation may become very strong. 
When the optimal point happens to be near the critical surface it is quite likely that the 
aformentioned oscillations take us to points which lie on the outside. This will lead, of course, 
to a failure of the method. In view of this problem we have devised a different strategy of 
computation which does enable us to locate the desired minimum no matter how close it is 
to the critical surface. The relevant procedure is analyzed in detail in Appendix A. 

Another issue which is of importance to our analysis is whether the experimental data 
lead to an optimal point which lies inside or outside the critical surface. It would be of great 
use if we could know beforehand, so that we do not try to reach an elusive "inside" point. 
To this end, we have developed (see Appendix B) a method through which we can determine 
whether the minimum lies inside our outside the critical surface without actually evaluating 
its exact location. We can also estimate what percentage of the region, as defined by the 
experimental errors of the data (in particle ratios), is inside and what outside the critical 
surface. 

A final matter of methodological concern pertains to the relevance of Bose/Fermi statistics 
effects and the impact they might have on our analysis. As it turns out (see next section) the 
source of greatest worry, as far as discrepancies between predicted and experimental values 
are concerned, are the (negative) pion multiplicities, in particular for the S + S experiment. 
The problem is that the direct inclusion of quantum particle statistics is quite difficult to 
be accomodated by the SBM scheme [16], whose formulation is based on the Boltzmann 
approximation. So, in order to gain an approximate evaluation of the error that is due to 
the omission of the Bose/Fermi statistics we turn to the Ideal Hadron Gas (IHG), always 
with the use of the grand canonical ensemble. 

■^Recall that the SBM equations have analytic solutions only in the region inside the critical surface. 
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The method that we shall adopt consists of finding the optimised set of parameters 
(VT^ / An'^ ,T, \q, \s,^s) from the S-SBM theory for the S + S data, subject, of course, to 
the constraint < S >s-sbm= 0. Then we determine the set {VT^ /An^ ,T, Xq, X'g,'~fs) which 
represents an optimal point in the Ideal Hadron Gas (IHG), formulated in the Boltzmann 
approximation and in the grand canonical ensemble (a constraint of the form < S >ihg— 
is naturally imposed). For the latter collection of parameters we calculate the particle mul- 
tiplicities {NiHG)i that we are interested in. At this point we switch to Bose/Fermi statistics 
as per in the IHG. For this case (BF) and for the set of parameters (VT'^/47r'^, T, A^, A^',7s) 
which verify the constraint < S >bf— 0, we calculate, once again, the desired particle mul- 
tiplicities {NBF)i- From the last two evaluations it follows that the relevant error in the 
estimation of the ith multiplicity, when we neglect the quantum statistics in the IHG, is 
fi = [{NiHo)i ~ {NBF)i]/{NjHG)i- OuT working assumption is that about the same amount 
of error persists in the case of the SBM. So any correction factor for this situation will be 
taken to be ~ 1 -|- /j. This correction factor will only be used for the h~ multiplicity in the 
S + S experiment and, as it will be seen, is small (of the order of 3%). 



4. Analysis of the Experimental Data 

In this Section we shall present the results of the application of the methodology we just 
described to experimental data. Our primary emphasis will be placed on particle multiplici- 
ties extracted from the NA35 S + S experiment at CERN. We shall, for purposes of checking 
the consistency of our scheme, also confront data from the p+p experiment UA5 at CERN. 

We shall display our results on multiplicities in a series of Tables where a fit with 
respect to the optimal set of thermodynamical parameters will be given as well. A series 
of the graphical presentations of our results will be depicted on the (T, /ig) plane. We shall 
first set 7s to its value determined from the fit- We also drop the dependence on the 
volume V, otherwise we would have to make a specific choice before being able to display 
our results on a two dimensional plane. Making our plots as independent as possible on 
fitted parameters, we shall turn to particle ratios. The ratios that will be used will be chosen 
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via a procedure that will be described in Appendix C. For a given particle ratio x, with 
experimental value x^^^, we shall plot, for fixed 7^, the projections on the (T,fiq) plane of 
the fines x{T, fig, jis, Is) = x'^'^^ + 5x^^^ and x{T, fig, fig, 7^) = x*^^^ — Sx'^'^^. This will be done 
for each the ratio. On the same figure we shall trace the projection on the (T, /ig) plane of the 
intersection of the surface < 5" >= with the critical surface (calculated for the same value of 
7s), which sets the limit of the hadronic domain. Finally, the fitted values of T and /ig, which 
result from the multiplicities, will also be indicated, along with the errors that result from 
the fit, represented by a solid circle on each graph. In this way we can conclude whether 
all the bands that are formed from the ratios and their errors have a common overlapping 
region. Such an occurrence would verify that the measured multiplicities are consistent with 
a system in thermal and (partial) chemical equilibrium described by the S-SBM and would 
help us evaluate the corresponding T and fig. More interestingly, we can determine how close 
to the limit of the hadronic domain the common region of the particle ratios' bands is and 
to what extent that region lies either in part, or as a whole, outside this boundary. Finally, 
we shall be in position to compare the results from the fit and the graphical analysis from 
the particle ratios, so that we can evaluate the degree of their complementarity. 

In the case of p + p experiment similar graphs will also be given, with the difference that 
they will be depicted on the (T, 7^) plane, while no choice for one of the relevant parameters 
will be needed in order to draw the bands which correspond to the particle ratios. 

4.1. S + S experiment (NA35) 

The study of S + S experimental data, which pertain to an isospin symmetric nucleus- 
nucleus colliding system, will be conducted first in the full phase space and subsequently at 
midrapidity. 

4.1.1. Full phase space data 

We deal, first, with the multiplicities that cover the full phase space. The An data from 
the S + S experiment at SPS (200 GeV/nucleon) [23] are listed in Table 1, along with 
experimental errors. With B — B we denote the baryon — antibaryon multiplicity. In the 
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present attempt we confront the full set of data which includes h" (mostly tt") multiplicities. 
A separate analysis which excludes the /i~'s will also be conducted for reasons that will 
become obvious in the following. 

4. 1.1. a. Analysis with the h~ 

In this analysis, the free phenomenological parameter of the model Hq is fixed so that 
To attains, for 7^ = 1, its "maximum" value Tq — 183 MeV (see [19]). The theoretically 
calculated multiphcities are hsted in the third column of Table 1. For the theoretical estimate 
of the h~ multiplicity, we have used a correction factor / = 1.0255 to account for the effect 
of quantum statistics. The evaluation of this factor has already been elucidated in Section 
3. 

The values of the optimised thermal parameters {VT^/An^, T, A^, 7s) which are calculated 
from our fit are hsted in the second column of Table 3, along with their errors. Here we 
have to note that in the collision of two S nuclei it is not evident how many nucleons really 
interact, even for central collisions, so the baryon number is not known a priori. This can be 
seen from the B — B multiplicity in Table 1, which is the net baryon number per interaction. 
It is measured to be 54 ± 3, whereas the total baryon number of the two incident S nuclei is 
64. So jJLq has to be left as a free parameter to be determined by the fit. The aforementined 
set comprises the true collection of free parameters which remain on account of the < 5" >= 
constraint. Accordingly, the degrees of freedom (dof) of the fit are 9 — 4 = 5, given that 
we have 9 experimental points. The least value of which is achieved by this fit is also 
listed in the same column. For completeness we have included the value of A^, which does 
not belong to the independent set of parameters, as well as those of the chemical potentials 

and iig-i which result from the calculated fugacities. 

The analysis with the particle ratios is depicted in Fig. 5, where we have set 7^ = 0.664. 
The chosen particle ratios are listed in Table 2, along with their experimental values and 
errors; they will be employed throughout our treatment of the full phase space data for 
the S + S experiment. In Fig. 7 we depict the value of {Ng^p — Ntheory) / (^exp, that is the 
number of standard deviations between experimental and theoretical values. The particular 
way of displaying our results was so chosen because there is great difference in the order 
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of magnitude of the measured multiplicities. The points of our fit are denoted by empty 
squares. 



S + S (NA35) 47r phcise space 



Particles 


Experimental Data 


Calculated with h 


Calculated without h 
(Case B) 


K+ 


12.5 ±0.4 


12.581 


12.706 


K- 


6.9 ± 0.4 


7.4590 


6.6332 


Ks' 


10.5 ± 1.7 


9.8106 


9.3791 


A 


9.4 ± 1.0 


7.8385 


9.6767 


A 


2.2 ±0.4 


1.3720 


2.0156 


P 


1.15 ±0.40 


1.9994 


1.5117 


p-p 


21.2 ± 1.3 


22.849 


21.529 


B-B 


54 ±3 


53.544 


52.348 


h- 


98 ±3 


94.086^ 


71.227'' 



A correction factor 1.0255 has been included for the effect of Bose statistics. 

A correction factor 1.0171 has been included for the effect of Bose statistics. This multiplicity is not 
included in the fit. 

Table 1. Experimentally measured full phase space multiphcities in the NA35 S + S 
experiment and their theoretically fitted values by S-SBM, with the inclusion of the h' 

multiplicity and without it (case B). 

It is evident that the fit, where all the multiplicities are included, is not so good, as can 
be witnessed from the relatively large value of (16.73, with 5 degrees of freedom). This is 
similar to the result of Becattini {x^/dof = 17.2/5) [9] and of Sollfrank (x^/dof = 11.6/4) 
[10], who have conducted a fit to the multiplicity data from NA35 S + S (with h~ included) 
for the ideal hadron gas case, formulated in the canonical ensemble. 

Looking at Fig. 5 one can see that there is no common overlapping region between the 
particle ratios. The fitted value for the pair (T, Hq) seems to lie inside the hadronic phase, 
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albeit near its limits. An analysis conducted along the lines described in Appendix B leads, 
for the values of particle ratios listed in Table 2, to a probability of 92.97% (238 points out of 
256) for lying insideQ and 7.03% (18 points out of 256) for lying outside the hadronic phase. 
Another point of note is that the strange sector is greatly suppressed, leading to 7^ = 0.66. 



S + S (NA35) 477 phase space 



Particle Ratios used 


Experimental Values 


Calculated without h 
(Case A) 


K-/K+ 


0.552 ±0.037 


0.52811 




0.84 ±0.14 


0.74268 


A/K+ 


0.752 ±0.084 


0.73901 


A/K+ 


0.176 ±0.032 


0.14413 


p/K+ 


0.092 ±0.032 


0.12044 


p — pI 


1.70 ±0.12 


1.7191 


B - B/K+ 


4.32 ±0.28 


4.1593 


h-/K+ 


7.84 ±0.35^= 


5.8036'^ 



^ This ratio is not included in the fit. 
A correction factor 1.0199 has been included for the effect of Bose statistics. 

Table 2. Particle ratios from the experimentally measured full phase space multiplicities in 
the S + S experiment and their theoretically fitted values by S-SBM, without the inclusion 

of the h~ / particle ratio (case A). 

4.1.1.b. Analysis with the h excluded 

We now perform the same analysis, but without explicit input from the h~ multiplicity. 
The latter can, of course, be calculated theoretically but with the optimal parameter set 
which results from the fit to the rest of the multiplicities. It is now determined that, for the 
same choice of Hq, the values of the thermal parameters lead outside the hadronic domain. 
Not being able to locate the point where attains its absolute minimum in this case, we 
^This probabihty will be refered to as Pinside and its value will be hsted to the appropriate Tables. 
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search for that point on the critical surface which optimizes the value of . Clearly, this 
is the closest we can get to the absolute minimum. In order to carry out such an analysis 
we have to turn to the method described in Appendix B. Wc designate this procedure as 
case A. The theoretically calculated particle ratios are listed in the last column of Table 2 
and are plotted with empty triangles in Fig. 8. The parameters that represent the location 
of the aforementioned least value of are listed in the second column of Table 3. Values 
for the volume have not been entered in this column, since particle ratios were used in the 
analysis. Note also that the fitted parameters arc not accompanied by errors, which cannot 
be evaluated due to the fact that we are not at the absolute minimum value of x^. In Fig. 
6 we have plotted the bands for the particle ratios that have been employed in the analysis, 
while the single point designates the location of the least value of ^ in the hadronic phase 
(on the critical surface) . 



S + S (NA35) 47r pheise space 



Fitted Parameters 


Fitted with h 


Fitted without h 
(Case A) 


Fitted without h 
(Case B) 


T (MeV) 


169.7 ±9.1 


176.7 


184.3 ±5.1 


\ 


1.538 ±0.064 


1.609 


1.613 ±0.079 


Is 


0.664 ± 0.056 


0.909 


0.96 ±0.16 




1.61 ±0.51 




0.61 ±0.31 


xVdof 


16.73 / 5 


3.06/4'= 


2.62 / 4 


K 


1.089 ±0.056 


1.010 


0.989 ± 0.052 


l^q (MeV) 


73.1 ± 8.0 


84.0 


88.2 ± 9.3 


Pis (MeV) 


14.4 ±8.8 


1.7 


-1.9 ±9.7 


PiNSIDE 


92.97% 


25.78% 





*= It is the minimum of within the Hadron Gas with Tq = 183 MeV (for 7^ = 1), not the absolute 
minimum. 

Table 3. Results of the analysis by S-SBM of the experimental data from S + S experiment 
{Att phase space) , with the inclusion of the h~ multiplicity and without it (cases A and B) . 
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What can be seen from the last analysis is that the quality of the fit is very good. The 
value of (though not its absolute minimum) is very small (3.06 with 4 degrees of freedom), 
compared to the previous one. All the calculated particle ratios that enter the fit are in very 
good agreement with the experimental values, as can be seen from the second column of 
Table 2 and Fig. 8. The last ratio, however, which includes the h~ multiplicity is very 
far from its experimental measurement. One can draw the conclusion that, for the thermal 
parameters coming out of this fit, the theoretically calculated number for the /i~'s (mostly 
7r~'s) is much less than the experimentally recorded one. 

From Fig. 6 one can see that bands for all the particle ratios, except for the one with 
the h~ multiplicity, converge to the optimal point. If the hadronic phase were allowed to 
occupy more space, one could infer that a good overlapping region would be formed. We 
evaluated the probability of being outside and inside the critical surface. It turns out that 
there is a 25.78% chance for being inside (33 points out of 128) and 74.22% for being outside 
(95 points out of 128). Finally, one notices that the strange sector is almost fully saturated, 
as 7, = 0.909. 

To complete the analysis without the /i~'s we perform another kind of fit, referred to 
as case B. This time we allow for the hadronic phase to occupy more space, so that the 
point of absolute minimum falls inside this phase. Such an arrangement can be achieved 
via the unphysical situation wherein Hq is set so that Tq = 193 MeV^ for 7^ = 1. From 
the relevant study we can gain a better feeling of the quality of the fit without h~ and to 
what extend full strangeness is attained. Our results indicate that when we vary Hq, which 
amounts to shifting the critical surface, the minimum value of remains basically constant. 
All the fitted thermal parameters, except for T and VT^/An^, undergo very slight changes. 
Our numerical results for the fit of case B give /dof = 2.62/4 and 7^ = 0.96 ± 0.16, i.e. 
the strange sector is fully saturated. The calculated values of the multiplicities are listed in 
Table 1 and plotted as empty triangles in Fig. 7. One notices again the difference between 
the experimental and the theoretical values of h~. The optimised parameters for this fit are 
listed in Table 3 (last column). 



^With this value the absolute minimum is just enclosed in the hadronic phase. 
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4.1.2. Midrapidity data 

We now turn to the analysis of multiplicity data restricted to the midrapidity region. The 
data are taken from Ref. [23,3] and are listed in Table 4. The rapidity interval is 2 < y < 3 
for all multiplicities, except for the p which are measured in the interval 3 < y < 4. This 
poses no problem since the system is symmetric around y — 2>. 



S + S (NA35) midrapidity 



Particles 


Experimental Data 


Calculated with h 


Calculated without h 


K+ 


3.2 ±0.5 


3.7659 


3.2874 


K- 


2.2 ±0.5 


2.4957 


2.0594 


A 


2.05 ±0.2 


1.9184 


2.0593 


A 


0.57 ±0.2 


0.33815 


0.53988 


P 


0.4 ±0.1 


0.47255 


0.41052 


p-p 


3.2 ±1.0 


4.3518 


3.0375 


h- 


26 ± 1 


25.544-'^ 


16.617» 



f A correction factor 1.031 has been included for the effect of Bose statistics. 



^ A correction factor 1.021 has been included for the effect of Bose statistics. This multiplicity is not 
used in the fit. 

Table 4. Experimentally measured midrapidity multiplicities in the S + S experiment and 
their theoretically fitted values by S-SBM, with and without the inclusion of the h~ 

multiplicity. 

4.1.2.a. Analysis with the h~ included 

Once again we perform a fit with the h~ included. The theoretical multiplicities are 
listed in Table 4 and plotted in Fig. 11. For the h" a correction factor / = 1.031 has been 
used to account for the effect of the Bose/Fermi statistics. The optimal set of values for the 
thermal parameters of this fit are listed in Table 6. The quality of the fit is, once again, not 
so good: x^/dof — 5.59/3. This also becomes evident from Fig. 9, where we have plotted the 
ratios. From this graph we can see that there is no common overlapping area of the particle 
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ratios. The location of the minimum of is observed to be well inside the hadronic domain. 
We actually find 100% probabihty for the point representing the experiment to be inside the 
hadronic sector (64 points are inside out of 64). Strangeness is also found suppressed, as 
7, = 0.78 ±0.12. 



S + S (NA35) midrapidity 



Particle Ratios used 
for the fit with h~ 


Experimental 
Values 


Particle Ratios used 
for the fit without h~ 


Experimental 
Values 


K+/h- 


0.123 ±0.020 


K+/A 


1.56 ±0.29 


K-/h- 


0.0846 ±0.020 


K-/A 


1.07 ±0.27 


A/h- 


0.0788 ± 0.0083 


A/A 


0.28 ±0.10 


A/h- 


0.0219 ± 0.0077 


p/A 


0.195 ±0.052 


p/h- 


0.0154 ±0.0039 


p-p/A 


1.56 ±0.51 


P - p/h' 


0.123 ±0.039 


h-/A 


12.7± 1.3» 



s This ratio is not included in the fit. 



Table 5. Particle ratios from the experimentally measured midrapidity multiplicities in the 
S + S experiment, used in the analysis with and without the h~ multiplicity. 

4.1.2.b. Analysis with the h excluded 

A different situation from the one described above arises when the h~ arc excluded. The 
results from the fit are listed in Tables 4 and 6, whilst the particle ratios that have been used 
are exhibited in Table 5. All the multiplicities, except h~, are consistent with a system in 
thermal and chemical equihbrium as described by the S-SBM. This occurence can be directly 
inferred from the very small value of x^/dof = 0.172/2. 

In Fig. 10 we depict the situation for the particle ratios in the (T, Hg) plane. All the 
particle ratios, except the one that includes the h~ multiplicity, form a common overlapping 
area which is partly within the hadronic domain. We find that the probability for a point 
belonging to this area to lie inside the hadron phase is 81.25% (26 points are inside out 
of 32) and outside 18.75% (6 points are outside out of 32). Note, also, that strangeness is 
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almost oversaturated: 7s = 1.15 ± 0.26. Of course, the error is such that the calculated 7^ 
is consistent with the maximum possible value of 7^ = 1. Note, on the other hand, that 
SoUfrank et al. in [3] have found 7^ = 1.19 for the S + S data at midrapidity. The study 
used the ideal hadron gas formulated in the grand canonical ensemble, with the cut to the 
input mass spectrum set at 2000 MeV. In our study, we have used as input all hadrons with 
mass up to 2400 MeV. 



S + S (NA35) midrapidity 



Fitted Paraiiieters 


Fitted with fr 


Fitted AA'itliout // 


T (MeV) 


156.7 ±8.1 


171.7 ±9.4 


\ 


1.485 ±0.081 


1.450 ± 0.075 


Is 


0.78 ±0.12 


1.15 ±0.26 




0.62 ±0.15 


0.22 ±0.13 


xVdof 


5.59 / 3 


0.172 / 2 


K 


1.149 ±0.050 


1.045 ±0.055 


l^g (MeV) 


61.9 ±9.2 


63.8 ±9.6 


pis (MeV) 


21.7 ±6.9 


7.6 ±9.0 


PiNSIDE 


100% 


81.25% 



Table 6. Results of the S-SBM analysis of the NA35 S + S experiment (midrapidity), with 

and without the inclusion of the h~ multiplicity. 

4.1.3. General observations on the S + S data analysis 

Some general remarks are in order. Let us start by commenting on the correction factors 
due to Bose/Fermi statistics. This correction is worth taking into account only for the h~ 
multiplicity. To demonstrate our claim, we consider the magnitude of the error of omitting 
the correct statistics by evaluating it for each multiplicity taken from the S — S experiment. 
Moreover, we choose that set of parameters which leads to the greatest correction factor. 
The results appear in Table 7. 
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Particles (i) 


("IHG)i ^ ' 


K+ 


0.510 


K~ 


0.578 


-*- ^ s 


0.527 


A 


0.228 


A 


-0.273 




-0.034 


P-P 


-0.092 


B-B 


-0.017 


h- 


3.043 



Table 7. Evaluation of the largest error in the analysis of the S + 5* data due to the 
omission of the Bose/Fermi statistics. In the evaluation reference to the Ideal Hadron Gas 

model is made. 

The analysis for both the full phase space and midrapidity give a fit which is not so good 
when the h~ are included, whilst all other multiplicities are very well fitted when the h~ 
are excluded. We conclude that the thermal parameters (temperature, chemical potentials 
and 7s) are more accurate when evaluated without the h" multiplicity. The measured 7r~ 
from the experiment are found to be, in both lot more compared to the theoretical 

predictions. The strangeness suppression factor 7^ is also near its maximum value 1 in the 
absence of the hr multiplicity. 

Next, we comment on the fact that the fitted temperature consistently drops when the h~ 
are included in the analysis. One may think that this should not happen, since the increase 
in temperature should lead to the increase of particle multiplicity. This is true, of course, 
provided the volume of the system is held fixed. In the fit we have performed the volume is 
not known; it simply enters as a free parameter. Thus what really is important is not the 
actual but the relative magnitude of the different multiplicities. When the h~ are increased 
much more than the other multiplicities, then the ratios of the latter over h" are decreased. 
Such a situation corresponds to lower temperature and probably greater volume (which, we 
repeat, was irrelevant to our analysis). 
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A final remark concerns the fact that at midrapidity the evaluated quarkchemical poten- 
tial fig is found to be somewhat less compared to full phase space case. This is expected, 
because the midrapidity region has smaller baryon number than the fragmentation region 
[23] (The full baryon number in the latter case is included in the 47r region). What is not 
so expected is that the fitted temperatures are lower at midrapidity than in Att rapidity 
regime. This can be explained by the difference that exists in the fitted values of 75. The 
temperature and 7^ are strongly correlated: increase in 7^ leads to decrease of temperature 
and vice versa (see Fig. 2). The fitted 7^ is found to be greater at midrapidity and that leads 
to a drop of the fitted temperature. Had we, on the other hand, calculated the temperature 
in full phase space and at midrapidity with 7^ fixed to the same value for both cases, then 
we would have found greater temperature for the midrapidity case. 

4.2. p + p experiment (UA5) 

We shall, now, concern ourselves with another isospin symmetric system, namely p + p 
collisions. An analysis of the multiplicity data of the relevant UA5 experiment has been 
performed by Becattini in [8], based on a canonical formulation of the ideal hadron gas. Our 
analysis will be performed through the grand canonical formulation of the S-SBM. The data 
we shall use are taken from [24,8] and are listed in Table 8, for the different center of mass 
energies. 

The examination of this system calls for a qualitatively different approach as compared 
to the S + S case. In the present situation, a p always interacts with a p. Thus the total 
baryon number is identically zero. We can thereby set the additional constraint on our 
system specified by < S >= 0. The < 5" >= constraint still exists, as well. These two 
equations can easily be solved analytically for the fugacities and their solution gives 

Aq = A, = 1 <^ //^ = = , (29) 

no matter what the value of 7^ is. So we are only left with 3 free parameters to evaluate 
from the fit and no constraint to fulfill. These parameters are {VT^/An^, T,^s)- In a similar 
manner as in the previous subsections, we perform a fit. In Table 8 we list the theoretically 
calculated multiplicities for the different energies and have plotted them in Fig. 15. In Table 

23 



8 we also list the ratios that we depict on the {T, 7^) plane in Figs. 12, 13 and 14. Our fitted 
parameters are tabulated in Table 9. 



p + p (UA5) 4n 


phase space 


Particles 


Experimental 


Calculated 


Particle 


Experimental 




Data 




Ratios used 


Values 


y/s = 200 GeV 


Charged 


21.4 ±0.4 


21.379 


Ks' / Charged 


0.0350 ± 0.0043 


Ks' 


0.75 ± 0.09 


0.77497 


n j Charged 


0.0350 ± 0.0047 


n 


0.75 ±0.10 


0.77828 


A / Charged 


0.0107 ±0.0028 


A 


0.23 ±0.06 


0.19934 


^ j Charged 


0.00070 ± 0.00070 




0.015 ±0.015 


0.01251 






= 546 GeV 


Charged 


29.4 ±0.3 


29.392 


K^ j Charged 


0.0381 ± 0.0027 


Ks' 


1.12 ±0.08 


1.1325 


A / Charged 


0.0090 ±0.0019 


A 


0.265 ±0.055 


0.29660 


^ j Charged 


0.00170 ± 0.00051 




0.050 ±0.015 


0.01984 






v^i = 900 GeV 


Charged 


35.6 ±0.9 


35.472 


K ^ j Charged 


0.0385 ± 0.0038 


Ks' 


1.37 ±0.13 


1.4288 


TL j Charged 


0.0281 ±0.0057 


n 


1.0 ±0.2 


1.1426 


A / Charged 


0.0107 ±0.0023 


A 


0.38 ±0.08 


0.32629 


^ j Charged 


0.00098 ± 0.00056 




0.035 ± 0.020 


0.02236 







Table 8. Experimentally measured in full phase space multiplicities and particle ratios from 
the p±p experiment and corresponding theoretical predictions by S-SBM. 

From Figs. 12-15 it can be inferred that a generally good fit exists. The minimum values 
of are similar to those found by Becattini in [8]. The overlapping region of the particle 
ratios is well inside the hadronic domain, as expected. The only discrepancy comes from the 
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S~ multiplicity for ^/s — 546 GeV which has an experimental value much greater than its 
theoretical prediction. This does not necessarily show something interesting. Finally, from 
the fitted parameters we can see that the temperature is almost constant with respect to the 
center of mass energy and that 7^ appears to be almost constant as well. 



p + p (UA5) 47r phase space 



Fitted Parameters 


= 200 GeV 


= 546 GeV 


^ = 900 GeV 


T (MeV) 


159.0 ±6.5 


159.3 ±6.3 


153.2 ±6.9 


Is 


0.422 ± 0.058 


0.454 ±0.027 


0.481 ± 0.054 




0.296 ±0.051 


0.398 ±0.064 


0.56 ±0.10 




0.448 / 2 


4.398 / 1 


1.583 / 2 



Table 9. Results of the S-SBM analysis of the data from p + p experiment. 



5. Assesments and Conclusions 

In this paper we have performed an S-SBM study of particle multiplicities, recorded in 
very high energy collisions in which isospin symmetry holds. On the theoretical side, we 
adopted the thermal description of the multiparticle system under partial chemical equilib- 
rium conditions. 

The relevant costruction involved five thermodynamical parameters subject to the con- 
straint < S >— 0. Our first concern was to determine an optimal set of values for these 
parameters through a fit to the experimental particle multiplicities. This formed the basis 
of our numerical computations which had dual aim. On one hand to assess the viability of 
our thermal description and, on the other, to locate the source of the multi-hadron system, 
on the phase diagramme, i.e. inside or outside the hadronic domain. 

The most striking observation, pervaiding our overall analysis, concerns the role of (neg- 
ative) pion multiplicities. A drastic change transpires, with respect to a fit of our results, 
between the cases with and without the inclusion of pions, irrespective of adjustments made 
for quantum particle statistics. Full chemical equilibration is consistently achieved when 
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pions are not included. Typical is the case of the An S + S data analysis which shows 
compatibility with a value for 7^ close to unity, i.e. with complete strangeness saturation, 
when the pions are not considered, whereas 7s drops to a considerably lower value when 
they are. The implications from analysis involving particle ratios are similar. Consistency 
with an equilibrium thermal behaviour is far better when experimental input from the h" is 
excluded. A similar situation arises in thermal analyses based on ideal hadron gas models, 
see e.g. [10]. 

One would like to understand the situation better, e.g. whether it is the freeze-out 
process responsible, or it is a signal that the location of the source of the hadrons is outside 
the hadronic domain. The latter case brings out the one special feature of the S-SBM, namely 
that it presents a critical surface beyond which the hadronic phase gives its place to a new 
one. Our studies involving particle ratios have established criteria for identifying the possible 
location of sources for the multiparticle hadronic state beyond this critical surface. The 
viability of such a possibility, assuming reasonable values for the bag constant (i?^/^ = 235 
MeV) and/or the maximum critical temperature Tq — 183 MeV, has come up in relation to 
our analysis of data from the S + S experiment NA35 at CERN [23]. In particular, the case 
of the An particle ratios (excluding the h~) indicates that about 3/4 of the particle-emitting 
state resides beyond the critical surface. This information, despite its semi-quantitative 
relevance, cannot be overlooked or ignored. We evaluate it as a strong sign of reaching and 
entering into the sought-for deconfined region beyond the hadron gas. 

Perhaps a more quantitative way by which to address the same issue is through entropy 
considerations. It was clear from our analysis that the n~ measured by the experiment are 
much more than we calculated theoretically. Sollfrank in [10] has pointed out that these 
particles may originate from a high entropy phase [25-27]. 

In order to shed some light into the notable discrepancy between the observed high h" 
multiplicity in the S + S data and the one predicted by S-SBM, we examine the entropy 
content, since it depends mostly on the multiplicity. We calculate this quantity for the 
thermodynamical parameters obtained from the fit to the An S + S data, Table 3. We remind 
that this fit produces a minimum point just beyond the assumed limit of the hadronic phase 
(To = 183 MeV, B^/^ = 235 MeV). 
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Entropy 


S-SBM 


341.73 


IHG 


264.68 


QGP (m, = 150 MeV) 


482.41 


QGP (m, = 300 MeV) 


464.36 


QGP (m, = 500 MeV) 


436.15 



Table 10. Calculated entropy for the fitted parameters of the 47r 5* + 5* data, case B. 

Table 10 presents the results of the respective calculations for the ideal Hadron Gas, S- 
SBM and QGP formahsms. We note that the S-SBM value is within (71 - 78)% of the QGP 
one for strange quark mass in the range 150 — 300 MeV and tends towards it. In contrast, 
the IHG value is within (55 — 57)% of the QGP value and about 77% of the S-SBM one. 
Noting that the predicted to the observed h' multiplicity is about 73% of the experimentally 
recorded value, we may infer that the S + S interaction at 200 GeV/nucleon may have just 
(initiated) a deconfinement phase transition and the h~ excess comes from an early stage of 
a deconfined quark matter state[|. 

In future work we intend to extend the S-SBM even more, by including the electric charge 
fugacity, Aq. This will enable us to take fully into account the difference in the number of 
the u and d quark of the initial colliding nuclei. In this way the analysis of a greater variety 
of heavy-ion experiments, including Ph + Ph collisions, will become possible. 
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Appendix A 

As mentioned in Section 3, in performing our analysis we have to find a set of parameters 
{VT^ /iiT^jT, \q, Xsjjs), in the context of the S-SBM, that corresponds to the minimum 
value of x^. In doing so, through a multidimensional generalisation of the Newton- Raphson 
method, we come across oscillations of points produced by the method during the evaluation 
procedure, irrespective of whether these lie inside the critical surface. We are, therefore, in 
need of another method for locating the minimum of that fully takes into account the 
restriction of the region where our functions arc analytically determined. Of course we have 
to locate the minimum under the condition that the constraint < >= is fulfilled (recall 
that our way of fulfilling this constraint was by the introduction of an additional parameter 
which played the role of a Lagrange Multipher) . 

Due to the fact that the critical surface is defined by the equation (/^(T, Ag, A^, 7^) = 
In 4 — 1, there is no particular value of one of the parameters through which we can tell when 
the critical surface is approached. A combination of values of all the parameters is needed. 
The solution to our problem calls for dropping one of the parameters in terms of which we 
have mapped the space enclosed by the critical surface and replacing it by (p. We choose to 
replace the temperature T. As far as the fulfilment of the < 5" >= constraint is concerned, 
we propose to try to find the minimum without ever leaving the < S >= surface. That 
means that one more parameter will no longer be considered as free, preferably Xg. Its value 
at every point will be determined by the rest of the parameters. 

The method can be summarised as follows. Prom the original set of parameters we choose 
to describe ^ ^ function of (VT^/in^, y?, A^, 7^) when multiplicities are the experimental 
input and of (99, Xq,"fs) when the input are the ratios. In the following we consider the case 
of the ratios. We start by an initial point which is inside the critical surface, denoted by 
(f/?'^, A^, 7°). Such a point is easy to specify: it is enough for it to satisfy </7° < ln4 — 1. In 
order to determine the value of this initial point we have to find the corresponding 

fugacity A^. For this reason we solve the system of the two equations: 



^(r,A;,A„7°) = 9'' 



,0 



(30) 



<5>(r,A;,A„7°) = o. 



(31) 
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The above system is solved numerically through a proper Newton-Raphson method. Let us 
stress that our designation of a starting point (T**"*"*, Af to the method is of paramount 
importance. If, again, we are led to points outside the critical surface then the present 
method will not solve our problem. To this end, it is very useful to observe that the critical 
surface around the value oi /ig — is almost perpendicular to the (T, /ig) plane (see Fig. 
9 in [19]). That means that whether we are outside or inside the critical surface will be 
determined by the rest of parameters and not by A5. So it is not important what starting 
point we shall give to A^, e.g. a value A**'"'* = 1 will not cause any trouble. The starting 
value for the temperature, on the other hand, needs more caution. This can be given by 
solving numerically the equation: 

(^(r^*-*,A;,Af-*,7°) = <^V (32) 

Once the desired T and A^ are specified can be subsequently evaluated. It is noted that 
the described routine has to be repeated every time the function is evaluated. 

The subsequent steps involve the performance of a number of one parameter minimisa- 
tions. A lot of routines exist for this task (see for example [28]). We consider the function 
X^(99, A^, 7^), which means that we hold fixed the parameters A^ = A^ and 7^ = 7^, and 
proceed with its minimisation. Suppose we find that its minimum is located at (p — (p^. 
Then we minimise x^(<^^, A^, 7^) with respect to A^. Let A^ be the location of this minimum. 
We repeat the routine for x^(</?^, A^, 7^), this time with respect to 7^, which gives us 7^ . Of 
course, the point (v?^, A^, 7^ ) will not give the location of the absolute minimum, so that the 
whole procedure will have to be repeated until is efficiently reduced. Note that, during 
our trials to find minima for one of the parameters A^ or 7^, we are not in danger of getting 
out of the critical surface because (p is held fixed to a value less than In 4 — 1. When we 
try to find a better (p, with Xg and 7^ fixed, we can get out of the critical surface only if 
the method of one parameter minimisation gives a value to (p greater than In 4 — 1. Usually, 
when this happens it means that the absolute minimum of is outside the critical surface. 
The method described in the following appendix will certify this fact. 

The method exhibited above is very slow compared to the generalised Newton-Raphson 
one, especially if the parameter VT^/An^ is added (in the multiphcity case). But it will 
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locate the desired minimum, even if it is near the critical surface, provided, of course, that it 
lies inside this surface. It can also be used to guide us near to the optimised parameters of 
the minimum, so the oscillation of the "test" points of Ncwton-Raphson is suppressed and 
we can reach the minimum with no danger of getting outside the region of analyticity. 



Appendix B 

In this appendix we shall present the method which allows us to conclude whether a given 
set of experimental data leads to thermal parameters which define a point inside or outside 
the critical surface. Because this method is rather slow we shall use as experimental input 
particle ratios and not multiplicities, reducing, in this way, by one our free parameters. 

The method consists of, firstly, locating the minimum value of on the intersection of 
the critical surface with the < S >= surface. This intersection is no longer a single curve, 
as in [18,19]. Given that the new parameter 7^ has been introduced into the S-SBM, the 
intersection is described by two variables which we choose to be {Xq,js)- For a given value 
of this pair the corresponding A5 cr can be found by numerically solving the generalisation of 
eq. (49) in [19]: 

^'^i.lJj cr; Xg er; cr) 

= . (33) 

Z^2- exp[G{y, Xgcr, Xs cr, Is cr)] 

Then the value of multiplicities prior to decays, which enter the theoretical value of the 
ratios, is given by (see eq. (21)): 



dz 



1 z-2 



dX., 

5 _ ^^{y, cr , Xs cr-i Ts 

By 



thermal 



^'^iVi ^q cri Xs cri'~1s cri 



iy^i Tcri Xq cri Xs cri Ts cr) 



VT: 



-fi 



dz 



, Aj, . . .) 



5 ^ViHi ^q en Xs en cr) 

^ ' By 



0. (34) 



2; = 2 - exp[G(y, en K en Is cr)] 

Note that V and Tcr need not to be known since they cancel in the ratios. 

Knowing how to evaluate on the intersection we perform, as in Appendix A, as many 
steps as required, each of which consists of two succesive one parameter minimisations. In 
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the end, we obtain a point — (A^,7^) which corresponds to the minimum value of 
on the aforementioned intersection. 

To proceed further we think that, if the absolute minimum of lies inside the critical 
surface, then by moving away from A™ towards the interior there should be some direction 
where decreases. On the contrary, if the absolute minimum lies on the outside, the 
minimum value of within the hadronic domain should lie on its limits. Any move from 
this location, in any direction inside the critical surface, should lead to the increase of x^- 
With these considerations in mind, we perform the following check. We form a grid of 9 
points: A^^ — (A™ + 5^,7^ + Sj), where 5i — — (5i,0,5i and 5j — —52,0,^2, with 5i and 62 
small numbers. The points A^^ are taken on a surface very close to the critical one, which 
is almost parallel to it. Such a surface can be determined by (f — InA — 1 — S3, with ^3 a 
small positive number. The points A^^ are located within the hadronic phase, very close to 
A™ but in different directions from it. We evaluate for all A"^^ . If one of the points A^^ is 
found to lead to value of less than the one it has at A"^, then the absolute minimum of 
is considered to be inside the critical surface. If, on the contrary, all the points A^^ possess 
greater values for x^ than the one at A^, then we surmise that the absolute minimum lies 
outside. 

With the method we just described we can evaluate the probability that an experimental 
point, with its experimental errors, is inside or outside the critical surface. Explicitly, suppose 
that we have n ratios given by the experiment. One of these ratios, e.g. Xi, takes the 
experimental value xl^^±Sxl^^ and so is somewhere in the interval {xf^^ — Sxl^^, xl^^ + Sxl^^). 
By the same token, a second ratio X2 can lie in the interval (x^^^ — Sx^^^, x^^^ + Sx^^^) and 
so on. By taking one value for every ratio from its corresponding allowed interval, we form 
an n-valued point which is within the error margins of our data. This point can be fed as 
experimental input to x^ and with the previous method we can find out whether it leads to 
an absolute minimum inside the hadron gas phase or not. We can, similarly, consider other 
points, forming different combinations among the allowed values of each ratio. The ratio of 
the number of points that are outside (inside)to the whole number of points considered give 
the probability that the data lead us outside (inside) the critical surface. 

In practice it takes a lot of time to process one point. So we only enter two possible values 
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of every ratio at the limits of its interval. Thus the points that we consider as experimental 
input in ^ are: 

with i\ — ±1,^2 = ±1, = ±1. This accounts to processing 2" points for every 

experiment that provides n particle ratios. 



Appendix C 

In principle, it is equivalent to use as experimental input in a fit a number N of 
measured particle multiplicities or a number oiN — \ independent particle ratios. There are 
numerous combinations oiN — \ ratios that can be formed from the same set of multiplicities. 
The way we proceed in practice is the following. After performing the fit with the 
multiplicities from the S ^ S data, we conduct a number of fits using different sets 
of ratios. We found a great deal of difference among these results as both the optimised 
parameters and the minimum value of exhibited puzzling variations. 

Before we analyse the situation let us recall the basic reasons for which it is needed to 
deal with particle ratios. Firstly, to locate the possible overlapping region of the bands of 
the particle ratios drawn on the (T, /i^) plane. Secondly, to evaluate the least value of x^ on 
the two-dimensional intersection of the critical surface with the < S >— surface when the 
absolute minimum lies outside the critical surface. Thirdly, to estimate the probability that 
the experimental measurements lead to an "inside" or an "outside" point. It follows that we 
cannot avoid facing the question which set of particle ratios to use in a way that it leads us 
to results as close as possible to those represented by the multiplicities. 

Translating the multiphcities to a set of ratios does not change anything as far as the 
centroid values are concerned. What is changed are the accompaning errors. The errors 
associated with the ratios depend on the way the multiplicities are coupled together to 
form them. To clarify this point, let us suppose that we choose the multiplicity Xm with 
the greatest relevant error (e^ = Sx^/xm — max) and couple it with all the rest of the 
multiplicities. Let us also suppose that among them exists a multiplicity with very small 
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relevant error (e^ = 5xn/xn)- This multiplicity will now enter in the ratio Xn/xm and the 



relevant error will be Cnm = y^n + ^m ^ ^n- When the multiplicity fit is performed, x„ 
plays an important role in the evaluation of the fitted parameters, since it has a very narrow 
interval for allowed values. The fit has to arrange the free parameters so that x^'^^"^^ 
comes close to Xn and the contribution to from this multiplicity is not too great|. On 
the contrary, the big relevant error to the ratio that includes Xn has the effect that the fit 
does not really take the smallness of e„ too much into account. Therefore, in order for the 
results of the fit with the ratios not to come out too distorted, we have to pay attention not 
to upset the connection of every multiplicity to its error. One simple way to do this is to 
pick that multiplicity with the smallest relevant error and form ratios by coupling the rest 
of multiplicities with it. 

We have checked in practice that the fit to the ratios chosen with the abovementioned 
logic really lead to the same results (fitted parameters and value of x^) as for the multiplicity 
fit. In the case of the An data from S* + S* we find that the and the h~ multiplicities have 
the smalest relevant error (~ 3%). In order to have the same set of multiplicities to both of 
our fits (with and without h~) we choose the and form their ratios with the rest. In the 
case of the midrapidity data from S + S the has the smallest relevant error (^ 3.8%), 
so it was chosen to form the ratios with the rest, in the case of the fit with the inclusion of 
h~ . When h~ is excluded from the fit, the smallest error lies with A (~ 9.8%) and thus this 
multiplicity was chosen. In the same way the charged multiplicity was chosen for all the fits 
to p + p data. 



The dependence of on the relevant error can be seen if we arrange its terms as 
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Figure Captions 

Figure 1 (a) Intersections of planes of constant s-quark chemical potential /is with the 

critical surface (p{T, jiq, jig, 7s) = ln4 — 1 for Tq — 183 MeV (for 7^ = 1). 

(b) Intersections of planes of constant q-quark chemical potential //^ with the critical surface 

ip{T, i^g, ^is, 7,) = ln4 - 1 for Tq = 180 MeV (for % = 1). 

Figure 2 Variation of the critical temperature, Tq, for zero chemical potentials /ig and /is 
with 7s. Tq is set to be 183 MeV when 7s = 1. 

Figure 3 Projection on the plane (T, /i^) of intersections of planes of constant q-quark 
fugacity \q {fJ^q/T = 0.4) with the surface < S >= for two values of 75. Tq is set to be 183 
MeV when 7s = 1. 

Figure 4 (a) Variation of the projection on the plane (//q,/is) of the intersection of the 
critical surface and the surface < 5" >= with 7^, for two different values of Tq (for 7^ = 1). 
(b) Variation of the projection on the plane (T, /Js) of the intersection of the critical surface 
and the surface < 5* >= with 7s, for two different values of Tq (for 7s = 1). 

Figure 5 Experimental Particle Ratios in the (T, Hg) plane for 5" + 5" experiment measured 
in 4:71 phase space. 7s is set to 0.664. The point and the cross are from the fit with the 
. The thick solid line represents the limits of the hadronic phase (HG) as they are set by 
S-SBM. The other line corresponding to the ratio KJ^ /K^ lies at the domain of negative jig, 
for T > 100 MeV. 

Figure 6 Experimental Particle Ratios in the (T, ^q) plane for S + S experiment measured 
in 47r phase space. 7s is set to 0.909. The point represented by the solid circle corresponds 
to the location of the least value within the hadron gas of x^, without the h~ . The thick 
solid hne represents the limits of the hadronic phase (HG) as they are set by S-SBM. The 
other line, pertaining to the ratio jK^ , hes in the domain of negative /z^, for T > 100 
MeV. 
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Figure 7 Comparison between the experimentally measured multiplicities in An phase space 
and the theoretically calculated values in the fit with h~ and without h~ (case B) for the 
S + S experiment. The difference is measured in units of the relevant experimental error. 

Figure 8 Comparison between the experimentally measured particle ratios in An phase 
space and the theoretically calculated values without h" (case A) for the S + S experiment. 

The difference is measured in units of the relevant experimental error. 

Figure 9 Experimental Particle Ratios in the (T, fiq) plane ior S + S experiment measured 
in midrapidity. 7^ is set to 0.78. The point and the cross are from the fit with the h~. The 
thick solid fine represents the limits of the hadronic phase (HG) as they are set by S-SBM. 

Figure 10 Experimental Particle Ratios in the (T, /ig) plane for S + S experiment measured 
in midrapidity. 7^ is set to 1.15. The point and the cross are from the fit without the 
h". The thick sohd line represents the limits of the hadronic phase (HG) as they are set by 
S-SBM. 

Figure 11 Comparison between the experimentally measured multiplicities in midrapidity 
and the theoretically calculated values in the fit with h~ and without h~ for the S + S 
experiment. The difference is measured in units of the relevant experimental error. 

Figure 12 Experimental Particle Ratios in the {T,js) plane for p + p experiment for ^/s — 
200 GeV measured in An phase space. The point represented by the solid circle corresponds 
to the fit- The thick sohd fine represents the hmits of the hadronic phase (HG), as they 
are set by S-SBM. 

Figure 13 Experimental Particle Ratios in the {T,js) plane for experiment for ^/s — 
546 GeV measured in An phase space. The point represented by the solid circle corresponds 
to the fit- The thick sohd fine represents the hmits of the hadronic phase (HG), as they 
are set by S-SBM. 
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Figure 14 Experimental Particle Ratios in the {T, 7^) plane for p + p experiment for ^/s — 
900 GeV measured in An phase space. The point represented by the solid circle corresponds 
to the fit- The thick solid line represents the limits of the hadronic phase (HG), as they 
are set by S-SBM. 

Figure 15 Comparison between the experimentally measured (solid circles) particle ratios 
in 47r phase space and the theoretically calculated values (empty squares) for the p + p 
experiment and for ^/s — 200 GeV, 546 GeV and 900 GeV. The difference is measured in 
units of the relevant experimental error. 
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